> 

i-H 



^ 






X 



INFRARED RENORMALONS AND POWER CORRECTIONS 
IN DEEP-INELASTIC SUM RULES* 



Xiangdong Ji 

Center for Theoretical Physics 

Laboratory for Nuclear Science 

and Department of Physics 

Massachusetts Institute of Technology 

Cambridge, Massachusetts 02139 



(MIT-CTP-2381 HEP-PH /9411312 Submitted to: Nucl. Phys. B November 1994) 
> ■ 
O 

\Q ■ Abstract 



Infrared renormalons and l/Q^ power corrections in deep-inelastic sum rules 
C\j ! are studied. The renormalization of operators with power divergence are dis- 

cussed. The higher-twist terms in the operator product expansion are shown 
to account for the residual soft contributions survived from the Kinoshita-Lee- 
. Nauenberg type of cancellation in Feynman diagrams. The presence of some 

Qv . degree of arbitrariness in the twist separation allows one to define the most 



convenient higher-twist operators suitable for a particular non-perturbative 
CL|| method. The discussion is focused on the Bjorken sum rule, for which the 

Q_i' 1/Q^ corrections are considered on a lattice. 
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I. INTRODUCTION 

Deep-inelastic sum rules — Bjorken sum rule [U, Gross-Llewellyn Smith sum rule pf, 
Ellis- Jaffe sum rule 0, and a whole zoo of similar sum rules ||^ — are fascinating things to 
study, both from theoretical and experimental points of view. These sum rules form one of 
the clearest tests of perturbative and non-perturbative Quantum Chromodynamics (QCD), 
and tests of our ability to solve the theory as well. Some of the sum rules were derived before 
QCD was advanted, in current algebras that the asymptotic QCD turns out to respect up to 
logarithms. In QCD, they can be derived using the technique of Wilson's operator product 
expansion (OPE) or its equivalent. The imprints of QCD in these sum rules are radiative 
corrections, which can be calculated perturbatively, and the higher-twist corrections that 
involve high- dimensional hadron matrix elements. This paper is intended to address the 
issue on the separation of these two effects and some related aspects in the OPE. 

In order to be specific, I will center my discussion on the Bjorken sum rule, which 
presently is under active experimental study 0. The sum rule is for the nucleon's spin- 
dependent structure function gi{x,Q^) and the QCD version reads 0, 

Ijgfix, Q') - g-{x, Q'))dx = ^^^(^.(g^))^?^ - ±c,{a,(^Q'))^ + ..., (1) 

where p and n label proton and neutron, respectively, and qa = 1.257 is the neutron /5-decay 
constant. Cj(as((5^) are the coefficient functions that can be calculated in perturbation 
theory as power series in as(Q^), 

C(a,) = Y. Cn<. (2) 

n=0 

The second term in the right-hand-side of Eq. (1) is the twist-four correction in which (O4) 
is the hadron matrix element, 

{PS\Oil\PS) = 25^(00, (3) 

with 

O^ = ugF^^'^^^u - dgF^'^j^d, (4) 

where \PS) is the nucleon state with momentum P'^ and polarization S^ {{PS\PS) = 
2P°(27r)'^(5^(0)), and u, d, and F^'^ are quark and gluon fields. Both C4 and (O4) depend 
on a factorization (renormalization) scale /x^, however, the dependence cancels out in the 
product. For simplicity, I have neglected the target mass corrections and the twist-three 
contribution 0, both of which are present at order 1/Q^, but inessential for the following 
discussion. The ellipsis denotes still higher-twist contributions that are suppressed by more 
powers of 1/Q^- 

It is commonly accepted that the first (twist-two) and the second (twist-four) terms as 
well as higher-order terms in 1/Q^ in Eq. (1) are separately gauge invariant and finite, 
and hence the twist-separation is well-defined. Indeed, the leading-order coefficient func- 
tion C2(«s(Q^) is often calculated in terms of the virtual Compton amplitude for a zero- 
momentum quark state in dimensional regularization |H| , in which perturbative higher- twist 



matrix elements vanish identically, and the non-perturbative higher-twist matrix elements 
are known to be calculable, for instance, in lattice QCD by subtracting the corresponding 



perturbative matrix elements |^, or in the QCD sum rule method |T^. It turns out, however, 
this traditional view on the twist-separation is at best formal, and one in practice is forced 
to provide more details on the individual contribution. Eventually, the twist expansion has 
the similar degree of arbitrariness as the factorization of the coefficient functions from the 
operator matrix elements in individual term of the OPE []TT|. This surprising feature of 



the twist-expansion in deep-inelastic sum rules was first recognized and studied, to the au- 
thor's knowledge, by A. Mueller [|T2|, although similar studies in vacuum current correlation 



functions had begun much early and the results are much more well-known p3HT6|. 

This paper grows out of attempts to understand the l/Q^ effects in the deep- inelastic 
sum rules |T^ and Mueller's paper on ambiguity of the higher-twist contribution |T^ . Since 



I believe most of the physicists working in deep-inelastic scattering are unaware of this issue, 
I intend to be pedagogical in my presentation. Therefore, the discussions in the paper are 
often mixed with results that are already known in the case of vacuum correlation functions. 
I hope the references are sufficiently clear so the reader has no confusion about what is new 
and what is not. 

In section II, I discuss the infrared (IR) renormalons in the coefficient functions of the 
Bjorken sum rule. I start with an empirical observation and a calculation of a bubble- 
chain diagram to convince the reader that the perturbation series in fact diverges. Then 
I state some general results known about the properties of the renormalon singularities. 
In section III, I consider the renormalization of operators with power-like divergence and 
explain why the multiplicative renormalization of log divergence cannot define the operators 
uniquely. The conclusion is that pure non-perturbative operators in general do not exist 
and higher-twist operators depend on handling of the perturbative contributions. In section 
IV, I derive a cut-off version of the OPE by separating the soft and hard contributions to 
Feynman diagrams. The result shows that the coefficient functions can be defined without IR 
renormalons, but depending on a momentum cut-off, and the higher-twist operators depend 
on the same cut-off due to perturbative power contributions. However, the sum of the twist 
expansion is cut-off- independent. In section V, I argue that there exists a large degree of 
freedom in choosing the higher-twist operators. According to this, I define a twist-four 
operator for the Bjorken sum rule suitable for lattice calculations. The corresponding twist- 
two coefficient function is subsequently evaluated. The final section contains the discussion 
on the practical significance of the 1/Q^ corrections in the Bjorken sum rule, and summary 
and conclusion of this paper. 

II. INFRARED RENORMALONS IN COEFFICIENT FUNCTIONS 

The IR renormalons, in the Bjorken sum rule, begin with the calculation of the coefficient 
function C2{ot{Q'^)) in dimensional regularization and minimal subtraction scheme (MS). 



According to text books [0 , the coefficient functions Ci are the coefficients of local operators 
in Wilson's OPE, 






\c2{asm)J^M - ^C,{as)0,f, + 



(5) 



where J5/3 = 'U7/375M — d'jp'j^d is the isovector axial current whose nucleon matrix element 
is qa- Other quantities are defined following Eq. (1). Since Eq. (^ is an operator identity, 
it is valid when sandwiched in any state. The canonical way of calculating 6*2(0:5 (Q^)) is 
to apply Eq. (^ to a single quark state with zero momentum. It is easy to see that the 
higher- twist matrix elements vanish in dimensional regularization, and 6*2(05) is just the 
"ratio" of the quark Compton amplitude to the quark matrix element of the axial current. 
In fact, this is exactly what has been down in Ref. [§], where 6*2(0^) was calculated up to 
three loops, 

^2(0,) = 1- ^^ ' -3.58 ^^ ' -20.21 ^^^ M +..., (6) 

for three fiavors. Here the renormalization scale has been chosen to be fi"^ = Q^, and so 
there are no logarithms of type In Q"^ / jj^ in c^. Notice that the coefficient c„ is growing 
rapidly with n . An estimate for the fourth-order term gives — 130(as((5^)/vr)^ [|T^. As a 



matter of fact, there are indications that c„ may grow like n! for large n with a fixed sign! 
To see the rapid growth of c„, consider the contribution from a n-loop bubble-chain dia- 
gram shown in Fig. la. Using dimensional regularization and modified minimal subtraction 
(MS), I find, 

c„+i(fig.la) = --f (£)'^Fn 2(1 - e)r(2 - e/2) ^ Q^— + 7)^ 



47r A-k' 



k=0 



/6r2(2 - e/2)r(e/2) V"'' r((n + 1 - A;)e/2)r(l - {n + 1 - k)e/2) 



r(4-e) J T{l + {n-k)e/2)T{3~{n + 2-k)e/2) 



(7) 



where Cp = 4/3 and /9q = 5 — 2n//3 represents the contribution to the leading-order P- 
function from the one- loop gluon self-energy {nf is the number of fiavors). The symbol Fn 
refers to taking finite part of the expression, which turns out to be non-trivial. With the 
help of the symbolic program MAPLE, I get the following result for the rational part, 

n Fn[...] 

20 2.5761918 x 10^^ 

21 5.4100006 X 10^° 

22 1.1902018 X 10^1 (8) 

The irrational part that involves Reiman (^-functions and vr is small compared with the 
rational part, at least at small orders. Eq. (P) indicates that the contribution from the 
bubble-chain diagram grows like ~ —5 x n\ at large orders. 

One might argue that the large coefficients from this particular diagram might cancel 
with those from other diagrams at the same order. However, indications are that the large n 
behavior of c„ for the sum of all n-loop diagrams is qualitatively similar to the bubble-chain 
diagram. To see that, consider the Borel-transformed series for C{as), 






C(6) = E -^h\ (9) 

n=0 ^■ 

Then C(a;s) can be constructed from C(6) through 

/•oo 

C(a,) = co+/ C(6)e-^/°^ci6. (10) 







Then a fixed-sign n! increase is refiected by the (renormalon) poles of C(6) on the positive 
real axis on the h plane. From an analysis of the infrared physics in QCD, 't Hooft argued the 
existence of these poles [^ . As I shall explain later in more detail, it is the infrared behavior 
of QCD that leads to the appearance of the renormalon poles, and hence the divergence of the 
perturbation series. The locations of the poles at 6„ = 4:Tm/ (3q {n = 1,2,...) were identified 
by Parisi through a study of the IR divergence of the corresponding super-renormalizable 
theory at 4 — e dimension [|^. A more accurate description of the analytical structure of 



the IR renormalons was obtained by Mueller, who showed that the QCD radiative effects 



change the poles into branch points in the complex b plane [0. The leading branch point 
at 6i = 4:71 /Po gives rise to a large n behavior of c„. 




Cn+i ~ n\n^ hr ' (11) 



where 7 is related to P2, the two- loop beta function, and 71, the one- loop anomalous dimen- 
sion of the twist-four operator |T^ . 

The existence of the IR renormalons destroys the Borel summability of the perturbation 
series since the integral in Eq.(pilD is now ill-defined. Unlike Borel-summable series, the 
IR-renormalon-plagued series refiects the truly divergent nature of the expanding quantity 
in perturbation theory; the divergence has nothing to do with the choice of the expansion 
parameter. Since at large Q^ higher-twist terms are exponentially small compared with the 
terms in the leading- twist, and since different ways of regularizing the Borel integral in Eq. 
(p!0|) produce an uncertainty of order of 1/Q^, the the nature of the OPE clearly deserves a 
more careful examination |jl5|. 



III. COMPOSITE OPERATORS WITH POWER DIVERGENCE 

Composite operators in QCD have a special importance: The chiral symmetry breaking 
is characterized by a non-vanishing expectation value of ijjijj in the vacuum. The electro- 
weak couplings of quarks are through vector and axial-vector currents ip'^^ip and ip'^'^'j^il^. 
The trace anomaly of the energy- momentum tensor are proportional to F"^ = F'^^Fap-, and 
the axial anomaly is proportional to F'^^Fap. Finally, Wilson's operator product expansion 
employs myriad composite operators. 

Composite operators are badly divergent in field theory, and are not made finite through 
the standard renormalization of Green's functions [|T8|. Rather, they are renormalized indi- 
vidually and have their own renormalization scale dependence. For these reasons, composite 
operators are normally not physical observables. When they do appear as a part of a physical 
quantity, the renormalization scale is determined either by physics consideration or cancelled 
by the scale-dependence of other quantities. 



According to the divergence behaviors of composite operators in QCD, they can be 
classified into three categories (I disregard operators with vanishing physical matrix elements, 
which include BRST-exact operators and the equations of motion operators [pO[]): 

• Finite operators: These include vector and non-singlet axial- vector currents of quarks 
and their divergence like ijjmijj, the energy-momentum tensor, the angular momentum 
density, the topological charge density FF, etc. These operators are associated with 
the symmetries of the QCD lagrangian, and are physical observables of the system. 

• Operators with logarithmic divergence: All twist-two and twist-three operators have 
logarithmic, and only logarithmic, divergence. Here I restrict myself to just local op- 
erators. Let me remind the reader that the twist of an operator is defined by the 
difference of the dimension of the operator and its rank of spin in representations of 
the Lorentz group. Although the notion of twist arises from study of deep-inelastic 
scattering, no reference here is needed to this application. [All operators can be classi- 
fied in terms of twist because of the Lorentz symmetry.] The twist-two and twist-three 
operators have just logarithmic divergence because their matrix elements in physical 
states are dimensionless quantities. They cannot mix under renormalization with any 
operators of lower dimensions. However, they can mix among themselves, yielding 
a dimensionless mixing matrix. The renormalized operators depend on the genuine 
renormalization scale fi"^. 

• Operators with power divergence: Operators of twist-four and higher normally have 
power divergence. Exceptions include finite operators discussed above and operators 
with special symmetry properties. The simplest example of operators with power 
divergence is F"^. One might argue that since the energy- momentum tensor is a finite 
operator and F^ appears in its trace (trace anomaly), so {(3{g)/2g)F'^ must be a finite 
operator. This is true as long as one is talking about the difference of the matrix 
elements in the excited states of QCD and the vacuum. Since the vacuum energy- 
momentum density is not a physical observable, {(3{g)/2g)F'^ in the vacuum needs not 
to be finite (the proper normal ordering of the operator is always implied here.) 

Once an operator has quadratic or quartic or higher-order divergence, they can mix 
under renormalization with lower-twist operators. For instance, F'^ can mix with the trivial 
operator 1, and '?/;F'^'^ 7,^75'?/' can mix with ijj'y'^ilj, etc. As shall become clear soon, these 
mixings have important implications about the twist separation in the OPE. 

Let me consider the renormalization of power-divergent operators in dimensional regu- 
larization. Although one cannot do non-perturbative calculations in this regularization, the 



coefficient functions in the leading- twist are normally calculated in the scheme [^]. Since 
perturbative QCD does not have any mass scale, the mixing of higher-twist operators with 
lower-twist ones vanishes identically. [Integrals of type / d'^k/k"^ are taken to be zero.] Only 
logarithmic divergence appears in the matrix elements of the operators, which can be renor- 
malized in a standard way. Thus it seems that composite operators can be defined up to their 
logarithmic divergence and are devoid of any power-dependent perturbative contributions. 
This standard treatment of composite operators in dimensional regularization is in fact 
deceptive. In principle, infrared and ultraviolet physics are entirely different and shall not be 



mixed. However, the statement J d'^k/k'^ = 0, for instance, has a ridiculous imphcation that 
the ultraviolet divergence is cancelled by the integration in the infrared region! Not only 
the integration is not reliable in the low k"^ region due to non-perturbative physics, but also 
is finite. Thus the implied cancellation can never occur in reality. Actually, the ultraviolet 
power divergence should not be thrown away like this because they dictate the correct 
operator mixing. Furthermore, a perturbative calculation in the low k"^ region generates IR 
renormalons just like the coefficient functions discussed in the last section. Thus J d'^k/k'^ = 
is mathematically convenient, but physically unwanted. 

Indeed, in a study of the two-dimensional 0(N) non-linear sigma model, David found 
that the composite operators with power divergence are not well-defined in dimensional 
regularization flSl. He calculated the spin- wave condensate in the physical vacuum using 



the regularization (which is possible in 1/A^ expansion) and found that the result is not 
unique. His argument for the presence of the ambiguity in the composite operator is subtle 
and roughly goes like this: When e = 4 — d is kept finit in dimensional regularization, the 
higher-twist matrix element has a series of poles on the positive real-axis of the e plane 
starting from e = the degree of divergence of the operator. The poles are present because 
at these e, one cannot distinguish a higher-twist contribution from a lower-twist one. These 
poles make the limit of e — > ambiguous, depending on whether one takes Ime > or < 0. 

The interpretation of David's ambiguity is that composite operators in dimensional regu- 
larization have IR renormalons. They were not there originally, but appear when the scheme 
insists to treat the contribution in the soft momentum region in perturbation theory. Those 
perturbative soft contributions are needed to cancel the ultraviolet power divergence be- 
cause the relations like / d'^k/k"^ = are used. Consequently the power divergence is killed 
with a price: emergence of IR renormalon singularities. It is these singularities that cause 
the ambiguity in David's calculation. [An explicit example of the renormalons in matrix 
elements is provided in the next section in Eq. (20)] Thus in dimensional regularization, in 
addition to subtract the logarithmic divergence of an operator, one needs to specify how the 
IR renormalons are regulated. 

What if one uses other regularization schemes, such as the scheme with a momentum 
cut-off? In a cut-off scheme, composite operators are finite, as long as the cut-off is kept 
finite. The operators are also free of any IR renormalon singularities. Their physical matrix 
elements depend on two scales: the non-perturbative scale generated from the breaking of 
the scale symmetry in the non-perturbative QCD [^ and the perturbative scale that is 



served as a cut-off for both the log and power divergences. If the degree of divergence of 
an operator is do, its matrix element goes like the do power of the cut-off when the cut-off 
approaches infinity. In light of dimensional regularization, it is tempting to define a pure 
non-perturbative matrix element by subtracting off the corresponding perturbative matrix 
element. Although the cut-off dependence now disappears, the IR renormalon singularities 
are introduced through perturbative matrix elements. Thus a pure non-perturbative matrix 
element is again not well-defined. 

To summarize the above discussion, operators with power divergence cannot be specified 
uniquely after subtracting their logarithmic divergence in dimensional regularization. One 
must specify in addition how the IR renormalons are regularized, or, equivalently, how 
the limit e — > is taken. In a cut-off regularization, the operators are cut-off dependent, 
independent of how the logarithmic divergence is regulated. For example, j3{g)/2gF'^ has 



no logarithmic divergence at all but still is cut-off dependent. One important consequence 
of this is that there are no physically well-defined higher-twist contributions in the OPE. 

Since the physical quantity that is expanded in the OPE is well-defined, the ambiguity in 
defining higher-twist operators must be correlated with the infrared renormalon singularities 
in the coefficient functions. In the 2D non-linear sigma model, David showed that the 
different choices for regularizing the Borel integral and defining the composite operators give 
the same final result. Thus despite the existence of arbitrariness in defining the contributions 
of different twists, the spirit of OPE is unspoiled. 

In fact, the situation with the twist separation is analogous to the factorization of 
collinear singularities in the OPE [|T^. There the leading-twist contribution in the OPE 
can be written as a product of the coefficient functions and the hadron matrix elements. 
The coefficient functions depend on the factorization scale because they are obtained by 
subtracting off collinear divergence in a single parton scattering. The matrix elements are 
scale-dependent because the local operators are renormalized operators disposed of ultra- 
violet divergence. The factorization scale controls how much physics is considered to be 
perturbative and how much is non-perturbative. Their product is scale- independent. 

IV. OPERATOR PRODUCT EXPANSION IN A CUT-OFF SCHEME 

According to the previous two sections, the traditional recipe for the OPE is at best 
formal: The coefficient functions are intrinsically divergent and higher-twist operators are ill- 
defined after removing the logarithmic divergence. However, both problems have a common 
origin — the IR renormalons. The existence of the renormalons makes the notion of twist 
separation procedure-dependent, however, specifics cancel in the physical sum. Specifically, 
one needs to define the limit e — ;> in higher-twist operators consistent with handling 



of the Borel integral for the coefficient functions ||T5[. However, the approach is hardly 
useful in practical applications. First, no one knows how to calculate the non-perturbative 
matrix elements of higher-twist operators in dimensional regularization in QCD. Second, 
despite tremendous progress made recently in calculating the QCD perturbative series to 



higher orders [|T^], it remains a formidable task to calculate the exact behavior of a infrared 
renormalon series. Lastly, even all the above are possible, the regularization does not seem 
to have clear physical motivation. [The author noticed a number of other proposals in the 



literature for regularization the infrared renormalons |12,23 



To my opinion, a better approach to the OPE, which was advocated by the former 



ITEP group p^ , |25[| , and was first taken seriously by Mueller |]T6|, is to consider the twist 
expansion entirely from the point of view of Feynman diagrams. Imagine a set of Feynman 
diagrams completely define a physical cross section or a correlation function in QCD. Often 
in accompany with a high-energy subdiagram, there are soft subprocesses mediated by low- 
energy quarks and gluons, which result in soft (infrared) and collinear (mass) divergence 



in perturbative calculations [|1^. In many physical observables, the soft divergence cancel 



according to equivalent of the Kinoshita-Lee-Nauenburg theorem ||2^, and the collinear 
divergence can be factorized into hadron matrix elements, which in the end are replaced with 
proper non-perturbative ones [|TT|. The remaining part is calculable as a perturbation series 



in the strong coupling constant a^ and is called infrared safe. The classical examples of the 
perturbative-calculable quantities include the total rate for e^e~ — > hadron, the coefficients 

8 



in Wilson's operator product expansion, and radiative corrections to jet production, Drell- 
Yan cross section, etc. 

However, it was pointed out by Shifman, Vainshtein, and Zakharov (SVZ) that using 
"infrared-safe" as a criteria to decide a perturbative calculation is inadequate p4|. They 



argued, in the case of current correlation functions in the vacuum, that the part of Feynman 
diagrams involving low-virtuality quarks and gluons cannot be calculated perturbatively 
despite such a calculation contains no infrared divergence and the contribution might be 
small. A legitimate calculation for the part must be done non-perturbatively, by using, 
in their case, various vacuum condensates. These vacuum condensates correspond to the 
higher-twist terms in Wilson's OPE. Thus, the role of higher-twist contributions is to take 
into account properly the residual soft contributions in Feynman diagrams after the KLN 
type of cancellation. 

Thus to devise a twist expansion, one can start by separating the hard and soft contri- 
butions in Feynman diagrams. This is precisely what Mueller did for vacuum correlation 
functions |16[. In this way, he obtained coefficient functions that are entirely perturbative 



and free of IR renormalons, and the higher-twist operators that are free of ultraviolet diver- 
gence due to an explicit cut-off. Both contributions now depends on the twist-separation 
scale but the sum does not. Notice the separation scale here plays the same role as the 
prescriptions of regularizing the higher-twist operators in David's paper. In this section, I 
follow Mueller's approach in Ref. |]l6l and attempt to construct a twist expansion for the 



Bjorken sum rule. To motivate the approach, let me come back to the physical origin of the 
renormalon singularities in the perturbation series. 

The appearance of n! in the coefficient functions is due to loop integrations in small 
momentum regions in Feynman diagrams [|T6| , p7| . In fact, consider a one-loop Feynman 
diagram with the gluon momentum k (chosen to be the loop momentum). Replace the 
gluon coupling as(Q^) (again, /x^ = Q"^) with a running coupling as{k'^), 

a^(k^) = '^^ -^ (12) 

Then the small k'^ integration is questionable because where the perturbative expression 
for asik"^) breaks down. [Notice the existence of the familar Landau ghost pole at k"^ = 
Q^exp(— as((5^)/3o/(47r)).] On the other hand, when expanding as{k'^) in terms of as(Q^), 

asik^) = a^) Y.i-m'^^^^T ln"(^), (13) 

the n-th term represents a legitimate partial contribution from the n-th order Feynman 
diagrams. When it is inserted in the loop integral, the repeated partial integrations give 
n\. A simple examination reveals that most of the contribution comes from the region of 
the integration k"^ ~ Q^e"", where the logarithms become large. Therefore the divergence 
of the perturbation theory is caused by illicit perturbative calculations in the soft regions 
where perturbative quarks and gluons are ill-defined. Although order by order, perturbative 
contributions from these regions are finite, they diverge as n ^ oo irrespetive the value of 
Q^ . Thus, a physically motivated OPE shall incorporate the perturbative soft (infrared) 
contributions in the higher-twist terms. 



To identify such contributions in the Bjorken sum rule, consider the one-loop contribution 
to the coefficient function shown in Fig. 2. The one-loop integrations are free of ultraviolet 
divergence after the standard QCD renormalization. The external quark momentum is set 
to zero after the collinear singularity in the box diagram (Fig. 2d) is subtracted by the 



one-loop matrix element of the axial current [^. Now focus on contributions from the 
region where the gluon momentum is soft [k"^ is small). In this region, one can expand the 
integrand in k'^/Q'^ and keep the leading-order contribution. Applying an upper cut-off A^ 
on the gluon virtuality, I find the perturbative part of Ci from the soft region. 



^^soft pert. ^ 8g^^ f^'' _d^ 



3g2 Jo (27r)4 fc2 

■9 TT Q2' 



(14) 



where in the second line have rotated the integration to the Euclidean space. The factor 
8/3 in the ffist line is the sum of 2, 1, and —1/3, coming from diagrams a), b), and c) 
in Fig. 2, respectively. The cut-off A^ should be on the order of 0.5 to 1 GeV^, above 
which perturbative calculation is justifiable. Clearly, the contribution in Eq. (0) must 
subtracted from ci and the corresponding non-perturbative contribution should be added as 
a higher-twist contribution. 

To find the correct non-perturbative contribution, I re-examine the one loop diagrams 
in Fig. 2. When the gluons are soft, we cannot contract the field operators in the per- 
turbative vacuum, instead we must keep them and evaluate the matrix elements in the 
non-perturbative zero-momentum quark state. [The same thing can be said of any quark 
lines that carry soft momenta.] For example, the diagrams (b) and (c) in Fig. 2 contribute 
to the quark Compton amplitude, 

where k and ki are quark and gluon momenta, respectively, with Euclidean kf restricted to 
< A^. H^'^°' represents part of the diagrams above the dashed lines and is perturbative. Sq 
represents that below and can be obtained from the ffist-order calculation of 

S.{k,k^) = J d'^d%e''<e''-^^{qm^i^{OA.{^,)^pmq{0)), (16) 

where ip and A^^ are quark and gluon field operators and \q{0)) is a zero- momentum 
quark state. Thus the one-loop diagrams arise from a perturbative expansion of the zero- 
momentum-quark wave function, which of course is incorrect as the gluon momentum be- 
comes soft. However, so long as one refrains from such a perturbative expansion, one gets 
the correct contribution from the soft-gluon region. 

To proceed, I use the technique of the collinear expansion as explained in full details in 
Ref. [^. The result is. 



8qa 
9Q 

where. 



Tr = '2^^'-''''J^M^)\oif\<im (17) 
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Qcnt _ f d ki d k2 Mc Jic iki^i ik2-e2 

X m [lpl5^DA^,) ■ zZ^x(6) + ^e'''^^'p^ns-fp^D^{^^)^D,{^,)] ^(0) (18) 

where p and n are light-like vectors that are introduced as a part of the coordinate basis. 
Clearly, as A^ -^ oo, Olf' -^ O^p, the local twist-four operator appearing in Eq. (||). For a 
fixed A^, O'^^p is a non-local operator. It is simple to show that at the one- loop order the 
perturbative calculation of the twist-four matrix element reproduces the soft perturbative 
contribution to the coefficient function in Eq. ([T^). 

The above result tells us several things about the operator product expansion in Eq. 
(5). First, the higher-twist operator has an upper cut-off in the Fourier components for the 
gluon fields, and thus is finite. To fully include the non-perturbative effect, the cut-off A^ 
must be larger than any non-perturbative scale, such as Aqcd- Second, when quarks and 
gluons have virtuality greater than A^, their effects are perturbative and are included in the 
coefficient functions. The separation of the two effects are scale-dependent, but the sum 
is not. Lastly, although I discussed only one- loop diagrams with one soft-gluon, the non- 
perturbative twist-four contribution takes care of all Feynman diagrams with one soft-gluon. 
Diagrams with one soft-quark line do not contribute at the twist-four level p^. Diagrams 
with two soft-gluon or quark lines are related to the twist-six terms, which are beyond the 
scope of this paper. 

Clearly, when taking O'^^ as a definition of higher-twist operator, the corresponding 
coefficient function in Eq. (^) shall be modified, 

icutx 



8 , N (^4/3 )pcrt 

-C4(a.)^ 



C2(«s) = C2pert(«.) + -C4(a.) ' '%" (19) 



where C2pert(as) the usual perturbation series in Eq. (|), C^^i^as) is the coefficient function for 
the twist-four operator. The perturbative matrix element of OZ* is to be evaluated in a zero- 
momentum quark state and is expected to have IR renormalons. Since soft contributions 
are subtracted in Eq. (|T^), we expect that C2{ois) is free of the leading IR renormalon 
singularity. However, the renormalon poles at Aim/ (3q with n > 2 are still present. To 
subtract them, one must consider the two- loop contributions to C2pert(tts) and the twist-six 
operators. 

To see that the perturbative matrix element of O'^^ has IR renormalons which cancels 
the leading IR renormalon in C2pert(tts), I consider again the bubble-chain diagram shown 
in Fig. la. The corresponding contribution to the matrix element is shown in Fig. lb. A 
simple calculation shows. 



(o^r)p-t(Fig.ib) = '^^cf 




-e{n-k)/2 

2\ e{n-k)/2 

(20) 

Thus the subtraction is both A^ and Q^ dependent. If I fix A^ at 1 GeV^, the subtraction 
is maximum at Q^ = 1 GeV^. As Q^ increase, the subtraction becomes smaller, indicating 
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less contributions in Cn from the non-perturbative region. For large n, the subtraction is a 
exceedingly slow-varying function of Q^. In fact, defining y = A^/Q^, I find at n = 20, 



y 


2Fn[...]/(3g2) 


0.1 


2.5761951 X 10^^ 


0.01 


2.5761877 x 10^^ 


0.001 


2.5755872 x 10^^ 



(21) 

Comparing with the result from Eq. (8), the subtraction is accurate up to one part per 
million for y = 0.1. As n — > oo, the large perturbative contributions to C2pert are entirely 
subtracted, independent of Q^ and A^. 



V. TWIST-FOUR CONTRIBUTION IN POINT-SPLITTING REGULARIZATION 

The cut-off scheme discussed in the previous section best illustrates the goal of the OPE: 
separating the perturbative and soft (infrared) contributions in a process. In analogy with 
the factorization of the coUinear singularities at each twist, the twist separation is not unique. 
In the cut-off scheme, this is shown by dependence of higher-twist operators on the cut-off 
scale A. The size of A must be large enough so that all non-perturbative physics is covered 
by higher-twist matrix elements, and at the same time shall be small enough that the pure 
perturbative physics is included in the coefficient functions of the leading-twist terms. In 
QCD, one believes that such scale exists around 1 GeV^. 

In view of what twist-expansion accomplishes, I argue that there exists a large degree 
of arbitrariness in choosing a cut-off version of the twist-four operator 04^: The specific 
cut-off scheme that I discussed before is only one example. As long as it fully incorporates 
the non-perturbative physics, as long as it approaches 04^ when cut-off is let go to infinity, 
any finite, non-local operator with right dimension and quantum numbers is a good twist- 
four operator. Once a choice is made, the coefficient function of the twist-two term can be 
calculated through applying the OPE to, say, a perturbative zero-momentum quark state. 
Thus the computation of the coefficient functions goes hand in hand with the choices of the 
higher-twist operators. 

The freedom of twist-separation is very useful in practical applications. The twist-four 
operator in the cut-off scheme is quite complicated, not suitable for non-perturbative meth- 
ods. In what follows, I seek a definition of the twist-four operator convenient for lattice 
calculations, because the lattice is the only practical means at present for a non-perturbative 
calculation. Following the definition, I calculate the corresponding twist-two coefficient func- 
tion. The approach here can be generalized straightforwardly to the application of the OPE 
to the QCD sum rules, which will be discussed in a separation publication [pO| . 

I choose to define composite operators in Euclidean space, in which the lattice realization 
is usually made. Imagine the space-time is discretized into a lattice with hypercubes; the 
lattice spacing is vr/A. The quark fields are assigned at every lattice point and gluon fields 
are represented by gauge links with 

U^{n) = pe"'^(^/^) /o' ^M("+^A)rfx (22) 
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representing a link, starting at the discrete point n and ending at a neighbouring lattice 
point in the direction (x. V is the path-ordering operator. I define an operator on lattice, 

04a(A^^) = ^J-^ea|3^luiJ{n)[U^{n)Uu{n + fi) - Uy{n)U^{n + i))]-ipil){n + jl + z>) + h.c. (23) 

where h.c. represents the hermitian conjugation term. The definitions of various Euclidean 
variables follow those of Rothe [0. It is simple to show that in the limit of A^ ^ oo, the 



operator goes to 0^^- Under the standard lattice gauge transformation, 

■ip{n)^ u{n)il){n) 
U^{n)^ u{n)U^{n)u\n + fi) (24) 

where u{n) is an element of the color SU(3) group, the operator is gauge invariant. Since the 
gluon fields are integrated along a straight path in a distance of vr/A, the large momentum 
components of the fields have effectively been cut off — the operator is finite. 

A few clarification is in order. First of all, I have said nothing about lattice QCD so far. 
I use the lattice as a tool to motivate a non-local, gauge invariant twist-four operator. Thus 
the actual quark and gluon fields and the location n'^ of the plaquets can still be regarded 
as continuous. [The method of defining finite operators is very similar to the well-known 
Schwinger's method of splitting the fermion fields in a current.] Second, the operator is 
non-local in Euclidean time, thus one has to rotate all loop integrals to Euclidean space in a 
perturbative calculation. Lastly, the operator has no complete Lorentz symmetry. However, 
the violation occurs at the scale of A around which physics is perturbative. When its matrix 
elements are used in conjunction with the coefficient function defined below, the symmetry 
breaking effects largely cancel. 

With the above definition of the twist-four operator, one can calculate the non- 
perturbative matrix element in a polarized nucleon state in lattice QCD. For instance, by 
choosing a series of lattice spacings a„ = 7r/(nA) (n = 1, 2, ...), one can construct the opera- 
tor from lattice variables almost trivially. By comparing a series of calculations with different 
n, one can extract the matrix element in the continuum limit. [Note that the scale A in the 
operator shall not be changed as the lattice spacing a approaches zero. The presence of this 
scale ensures the existence of the finite matrix element in the limit.] The lattice calculation 
is currently been persued by the MIT lattice group WI . 



What remains is to calculate the corresponding coefficient function of the twist-two 
operator. As usual, I sandwich the OPE in a zero-momentum quark state and calculate 
everything using dimensional regularization, including the perturbative matrix element of 
the twist-four operator. At one-loop level, the matrix element is. 



^"-^mt'^t'w- (25) 



7r2 



The calculation is done with Feynman rules in Euclidean space. From this we get the 
one-loop contribution to the coefficient function. 
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^1 = -r(l - T^TT^) (26) 



1 16A^ 

^*^ 97rg2 

The power term here is quite close to the power term in the cut-off scheme in the previous 
section. I emphasize again that this coefficient function must be used together with the 
specific twist-four operator in Eq. (^31) . 

Another way to define the twist-four operator is to let the non-locality of the operator 
in Eq. (|23|) A be the same as the lattice spacing a, approaching zero in the continuum limit. 
Then the lattice matrix element diverges like 1/a?. A finite twist-four operator can be 
defined by subtracting from the original operator its perturbative matrix element on the 
lattice, with a low momentum cut-off in loop integrals. Thus the twist-four contribution is 
again free of the IR renormalons but still depends on the infrared cut-off. Unfortunately, it 
is difficult to do perturbative calculations with a nucleon state. However, the idea is clearly 
applicable to vacuum condensates in the QCD sum rules [§,0]- 

VI. COMMENTS AND CONCLUSIONS 

According to the discussions in the previous sections, the leading-twist contribution to 
the Bjorken sum rule can still be regarded as a power series in as{Q'^)- However, the 
coefficient in each order is now not just a numerical number; it contains the subtraction of 
the soft contributions that are represented by a power series in IS? jQ^ . The most important 
soft-subtractions are those nominally suppressed by one power of l/Q"^- As the order of the 
perturbation n -^ oo, the subtraction grows like n!, independent of the value of Q^ and 
the subtraction scale A^. Therefore, the series in as would converge if other renormalons 
singularities and sources of divergent factors did not exit. [It is known that ultraviolet 
renormalons and the singularities in the Borel plane caused by instanton-anti-instanton 
pairs still make the series asymtotic, but I ignore them here.] 

Thus to the accuracy of l/Q^, one can successively take into account the QCD corrections 
in the following way. First calculate the coefficient function (72(0;^, A^/Q^) in powers of a^. 
At each other, the 1/Q^ subtraction must be calculated explicitly. At lower orders, the 
subtraction is smaller than the pure loop contributions. However, the former must be kept 
to maintain the accuracy. As n increases, both the loop contributions and subtractions 
decrease initially, and their difference also decreases. If the difference becomes the size of 
the non-perturbative 1/Q^ terms, the latter must be included. With the further increase 
of n, both the perturbative contributions and their subtractions start to grow. However, 
the cancellation now becomes more complete and the residual remains small. The first few 
terms of the Bjorken sum rule is. 



4 (04^A2)) 



97rg2^ \tx ) ' 'g2 

2 +••• (27) 



27 Q 

The soft contribution in the two-loop diagrams, denoted by r/, needs to be calculated. 

The non-perturbative twist-four contributions can be calculated on lattice, as explained 
in the last section. In the previous studies, the twist-four matrix elements have been evalu- 
ated in the MIT bag model and the QCD sum rule method. The bag model calculation has 
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no explicit cut-off dependence 0, which is a common problem for model estimates. Empir- 
ically, one may regard the cut-off being somewhere from 0.5 to 0.8 GeV^ |33|. The QCD 



sum rule calculations for the higher- twist matrix elements are clearly questionable |jTO|]. As 
I explained in sec. Ill, the operators used in these calculations are not well-defined. 

The IR renormalons exist in all perturbation series in QCD. The issue is closed related to 
the non-perturbative power corrections. Unlike the case of QED, the power corrections are 
an important subject in perturbative QCD. For instance, a precision determination of the 
strong coupling constant clearly needs a correct understanding of the higher-twist effects. In 
fact, already at present, the effects become the dominating theoretical uncertainty in some 



of the data analysis [0]. The framework described in this paper provides one way get the 
problem under control. 

To summarize, we have studied in this paper the IR renormalons and l/Q^ corrections 
in deep-inelastic sum rules. The result shows that the non-perturbative l/Q^ corrections 
depend on specific choices of the regularized higher-twist operators. I have proposed to define 
a twist-four operator for the Bjorken sum rule, which is suitable for lattice calculations. 
The coefficient functions for the leading-terms can only be calculated after the higher-twist 
operators are defined. Although most of the discussions in the paper are made for the 
Bjorken sum rule, all the results are clearly applicable to all deep- inelastic sum rules. 

I wish to thank A. Mueller for many useful discussions, comments, and suggestions. Dis- 
cussions with K. Johnson, W. Lin, J. Negele, B. Schreiber, U. Wiese are also acknowledged. 
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FIGURES 

FIG. 1. a). A bubble-chain diagram for the coefficient function of the Bjorken sum rule. b). 

The corresponding matrix element of the twist-four operator. 

FIG. 2. The complete one-loop Feynman diagram for the the coefficient function of the Bjorken 

sum rule. 
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